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ABSTRACT 

We present the results of a search for X-ray variability in the galactic open cluster 
NGC 2516. This cluster has been observed on 8 occasions using XMM-Newton mak- 
ing it an excellent target for such a study. We find 4 sources which show significant 
variability implying the detection of one significantly variable source every 25ksec. 
All sources exhibit spectra which can be fitted using an absorbed one or two temper- 
ature thermal plasma model. The brightest of these flares also show a prominent Iron 
line near ~7.0keV. All 4 sources lie near the observed main sequence of NGC 2516 
suggesting they are cluster members. We propose that at least 3 of the 4 objects are 
RS CVn systems. We compare the properties of the brightest flare with those of solar 
flares. 

Key words: Sun: flares - stars: activity, coronae - Galaxy: open clusters and asso- 
ciations: individual: NGC 2516 - X-rays: stars, binaries 



1 INTRODUCTION 

NGC 2516 is a rich, well studied, galactic open cluster. It 
has a low metallicity (eg Cameron 1985, Jeffries, James & 
Thurston 1998), is nearby (~387pc; Jeffries, Thurston & Pye 
1997), is young (~ 10 s years; eg Feinstein, Marraco, Mirabel 
et al 1973) and has a relatively low line of sight extinction. 
Because of its low metallicity, it is an excellent target to test 
how coronal activity and stellar structure depend on metal 
content. 

NGC 2516 has been the subject of several studies in 
the X-ray band. It was observed using ROSAT by Dachs & 
Hummel (1996), Jeffries, Thurston & Pye (1997) and Micela 
et al (2000). More recently it was the subject of a study 
using XMM-Newton (Sciortino et al 2001) and also Chandra 
(Harnden et al 2001). These studies have concentrated on 
determining the cluster X-ray luminosity function: Sciortino 
et al (2001) found that dwarf G and K stars in NGC 2516 
were less luminous than those in the Pleiades (which has 
a similar age but solar metallicity ratio). In contrast, the 
luminosity function of dwarf M stars in NGC 2516 and the 
Pleiades were indistinguishable. 

In this paper, we show the results of a search for X- 
ray variable objects in NGC 2516 using XMM-Newton data. 
With the relatively large number of X-ray objects it is an 
excellent target to determine the frequency of stellar flaring 
activity in a young, metal poor cluster. Further, NGC 2516 
has been observed on 8 occasions with XMM-Newton since 
it has been used as a calibration target. 



Date 


Orbit 


Obsld 


Clean Time 
(ksec) 


2000-04-06 


060 


0113891001 


17.6 


2000-04-06 


060 


0113891101 


15.0 


2000-06-10 


092 


0126511201 


26.1 


2000-10-31 


164 


0113891201 


0.5 


2001-01-29 


209 


0134531201 


17.0 


2001-01-29 


209 


0134531301 


5.0 


2001-06-24 


272 


0134531301 


0.0 


2001-10-29 


346 


0134531501 


16.7 



Table 1. The observation log of XMM-Newton observations of 
NGC 2516. We show the XMM-Newton orbit number and the 
observation id. 



2 OBSERVATIONS 

The satellite XMM-Newton was launched in Dec 1999 by 
the European Space Agency. It has the largest effective area 
of any imaging X-ray satellite. The EPIC instruments have 
imaging detectors covering the energy range 0.15-10keV 
with moderate spectral resolution. The observation log is 
shown in Table 1. The first two observations were centered 
at 7 h 58 m 20 s , -60° 52' 13" (2000), and the rest at 7 h 58 m 
22", -60° 45' 36" (2000). (We do not discuss the data taken 
using the Reflection Grating Spectrometer since the sources 
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Figure 1. The background subtracted X-ray light curves for each of the 4 sources which showed significant variability. The bin size is 
120 sec and covers the energy range 0.2-10.0keV. 



were comparatively faint and source confusion would be a 
problem. The Optical Monitor was in blocked mode). 

The data were processed using the XMM-Newton Sci- 
ence Analysis Software (SAS) v5.3.3. An image was gen- 
erated from each dataset in the energy range 0.2-5.0keV. 
We removed data at the edge of the detector furthest from 
on-axis since these rows had increased noise. We also se- 
lected those events which were tagged as FLAG=0: this ex- 
cludes events close to, for instance, hot pixels and chip 
edges. Time intervals which showed high levels of particle 
background were excluded. The resulting images were then 
source searched using standard SAS tasks. Events were then 
extracted around a specified radius from the source posi- 
tions. The events for each source were then subject to a 
Discrete Fourier Transform in the time interval of lmin to 
0.7xExposure length. All power spectra which showed a 
peak above a certain criteria above the mean power were 
flagged and inspected by eye. We also had a second flagging 
procedure which indicated which objects had a count rate 
which varied over a certain value. 

After the completion of this process, 4 sources were 
found to show significant variability in their X-ray light 
curves. Their X-ray positions were compared to the list of 
X-ray sources in NGC 2516 and their corresponding opti- 
cal counterparts (Jeffries, Thurston & Pye 1997). We found 
counterparts for all 4 sources: the offset between the XMM- 
Newton and optical counterpart was typically better than 
6" (which is the FWHM of sources in the EPIC pn). We 
show their X-ray positions along with their counterpart of 
the X-ray catalogue of Jeffries et al (1997) in Table 2. For 
each source we examined the other datasets to determine if 
the source was also detected at other epochs. In each case 
the source was detected but at low count rates (typically 
0.01-0.02 ct/s). However, they showed no significant varia- 
tion in these other datasets. 

We show the (background subtracted) X-ray light 
curves for each source in Figure 1. Three of the sources show 



JTP 


RA 


Dec 


V 


(B-V) 


Obsld 


80 


07 57 58.5 


-60 56 53 


12.84 


0.63 


0113891101 


138 


07 59 20.1 


-60 34 44 


11.96 


0.51 


0126511201 


114 


07 58 43.3 


-60 55 24 


13.97 


0.87 


0134531501 


122 


07 58 50.8 


-60 38 35 


9.51 


0.08 


0134531501 



Table 2. The X-ray sources which showed significant X-ray 
variability. JTP refers to the X-ray source number of Jeffries, 
Thurston & Pye (1997), the X-ray position (2000) is the position 
from the XMM-Newton EPIC images and the optical counter- 
parts are taken from Jeffries et al. We also indicate the obser- 
vation number (cf Table 1) in which the source variation was 
detected. 



evidence for a rapid flare, although in one of those (JTP80) 
the observation ends shortly after the start of the flare. The 
other source (JTP122) shows a more gradual rise in bright- 
ness. The most dramatic source is JTP138 which shows a 
gradual rise in flux lasting ~12 min followed by a sharp in- 
crease in flux taking place in less than 2 min. After reaching 
maximum, it reaches half its maximum intensity after ~20 
mins. The duration of the flare lasts around 90 min. We show 
in Figure 2 the light curves in the l-3keV and 3-10keV en- 
ergy bands for JTP 138 in time bins of 500 sec together with 
its softness ratio. Although the hardness ratio curve shows 
some evidence for a softening over the duration of the flare, 
it has large errors. However, there is some evidence for a 
small increase in flux in the l-3keV light curve just before 
the sharp rise that is not seen in the 3-10keV light curve, 
although again this may just be due to the low signal to 
noise. 



3 SOURCE SPECTRA 

We extracted spectra from each of the 4 variable sources: 
only the time interval of the flare duration and single and 
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Figure 3. The EPIC pn spectra of source JTP138. The best fit 
- shown as a solid line - is an absorbed two temperature plasma 
model. 
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Figure 2. The flare seen in JTP138 in the l-3keV (middle panel) 
and 3-10keV (bottom panel) energy band and the 1— 3keV/3- 
lOkeV softness ratio (top panel). 



double events were used. We also extracted background 
spectra from a nearby source free area off the detector. We 
used the 'canned' response files appropriate to the location 
of the source on the detector and generated an auxiliary file 
using the SAS task arf gen. We fitted the spectra using the 
X-ray fitting package XSPEC (Dorman & Arnaud 2001). 

All spectra could be fitted using an absorbed one or 
two temperature plasma model (we used the MEKAL model 
in XSPEC). We show the spectrum of the brightest flare 
(JTPI38) in Figure 3: the fit is good with no evidence for 
significant residuals at particular energies. This source also 
showed evidence for an Fe line at 7.0±0.2keV (the other 
sources did not). We show the spectral fits in Table 3. 

All spectra show a moderate amount of absorption, typ- 
ically ~5x I0 20 cm~ 2 , and temperatures of a few keV. In 
determining their luminosity we assume that all 4 sources 
are members of NGC 2516 and its distance is 387 pc (Jef- 
fries, Thurston & Pye 1997). In the case of JTPI38, where 
we observe the total duration of the flare, the energy re- 
leased was 8 x I0 34 ergs. In the case of JTP80, where we 
observed only the start of the flare, the energy released was 
greater than I x I0 34 ergs. There was a similar energy release 
observed in JTPf 14 where most of the flare was observed. 



4 THE NATURE OF THE SOURCES 

We have identified each of the 4 X-ray variable sources with 
an optical counterpart (Table 2). Jeffries et al (1997) identify 
JTPf 22 as having an AOV spectral type: the other 3 sources 



have not been classified spectroscopically. Further, JTPI38 
and JTPf 22 were both classed as cluster members by King 
(1978) from a proper motion survey. To determine the like- 
lihood of the other two sources being cluster members we 
plotted the V, B — V, data of Sung et al (2002) and traced 
the position of the main sequence. We then over-plotted the 
optical colours of our X-ray variable sources (Figure 4). All 
4 sources lie close to the main sequence of NGC 25f6: wc 
conclude that each source is likely to be a cluster member. 

Stellar coronal activity is common in stars which have 
a convective layer and a relatively rapid rotation period. 
Enhanced levels of activity are often seen in close binary 
systems, for example the RS CVn systems which typically 
consist of a giant or sub-giant primary with a main-sequence 
or sub-giant companion. In many RS CVn systems the ro- 
tation of one of the components is synchronised with the 
orbital motion, leading to enhanced rotation and high lev- 
els of activity, seen as large spot coverage and flares in X- 
ray, optical and radio bands. Additionally X-ray activity has 
also been seen on early-type stars, such as Castor A and B. 
Here it is thought the flaring component may originate from 
low-mass, late-type companions (Giidel at al. 2001). X-ray 
flares have also been reported in pre-main sequence stars (eg 
Schmidt f994) and single solar type stars (eg Landini et al 
1986). 

The energy released in the X-ray flares we observed in 
NGC 2516 were ~ x 10 34 " 35 ergs and lasted around 1-2 hrs. 
This contrasts to the energy released on a flare seen on Algol 
(7x10 ergs) which lasted over a day (Ottmann & Schmitt 
1996). In dwarf M stars flares have been seen in the U band 
(eg Byrne, Doyle & Butler f 984) although they are relatively 
rare and last several 10's of minutes. In X-rays, dwarf M 
stars show flares of energies typically ~ ^q 3 °- 34 er g S w ith a 
duration less than 1 hour (eg Cheng & Pallavicini 1991). We 
consider that it is unlikely that the flares originate from pre- 
main sequence stars since NGC 2516 is not a star forming 
region. Solar type stars have also been found to show flares in 
various energy bands with energies ~ 

1Q 33-38 ergg) although 

they are very rare (eg Schaefer et al 2000). 

In contrast, flares have been seen on RS CVn systems 
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JTP 


N H 
xlO 20 cm- 2 


kTi 
(keV) 


fcT 2 
(keV) 


Observed flux 
ergs s -1 cm -2 


Bolometric luminosity 
10 30 ergs s- 1 10 34 ergs 


* 2 
(dof) 


80 
138 
114 
122 


7+8 
'-6 

°-°-4.8 


3 1+2-1 
^•■■■-0.6 

n s+ - 4 

i o +0 - 7 

i ' U -0.5 


c Q+3.1 
^•'-1.7 


3.l±°;jj x 10- 13 
1.05±0.14 x 10" 13 
1.64X10" 13 
3.0±0.6 x 10- 13 


5 .5t; : ? 

4.9i?; 8 3 

s.oi;- 2 


1.2±0.2 
8.0+?-* 
1 l +o:2 

1 ' 1 -0.3 

°- D -i.o 


0.64 (39) 
0.97 (64) 
1.27 (16) 
1.67 (33) 



Table 3. The spectral fits to the four X-ray varying sources. We have extracted events covering only flare time interval. When determining 
the luminosity of the flare we assume a distance of 387pc and a solar metallicity. 



gives 2.7xl0~ 4 variable X -ray sources per ksec per X-ray 
source in field 1 and 3.0 x 10 -4 in field 2. Although simplistic, 
if every source in the field was a potential X-ray variable 
source, then these rates would imply that each source would 
show a flare every ~40 days. 



(B-V) 



Figure 4. We plot the location of the optical counterparts of the 
X-ray flaring sources in the V, B — V plane. We also show the 
location of the main sequence in NGC 2516 taken from the data 
of Sung et al (2002). 



which share similar characteristics to the flares we have ob- 
served, van den Oord, Mewe & Brinkman (1988) estimated 
the energy released in a flare on the RS CVn system a 2 
CrB was 2.4 xlO 34 ergs and lasted ~ 90 mins. Further, they 
found clear evidence for a softening of the light curve after 
the sharp rise in flux. This is comparable to the RS CVn 
system II Peg which released > 5 x 10 34 ergs in a flare last- 
ing more than several hours (Doyle et al 1991, Doyle, van 
der Oord & Kellett 1992). Their X-ray spectra are well mod- 
elled using a thermal plasma of a few keV. Further, there 
are several instances where Fe Ka emission has been seen in 
stars with active stellar coronae (eg Tsuru et al 1989, Giidcl 
et al 2001). 

We suggest that JTP 138, JTP 80 and JTP 114 are good 
candidates for RS CVn systems. The fact that JTP 122 has 
been classified as AO V spectral type may indicate that it 
has a dwarf M star as a companion. Followup spectroscopy 
will be able to test these assertions. 



5 THE FREQUENCY OF FLARING ACTIVITY 

NGC 2516 has been observed using XMM-Newton for 97.4 
ksec (cf Table 1). We detected 4 significantly variable X-ray 
sources over this time implying the detection of one X-ray 
variable source every ~25ksec. The first two observations 
were offset from the other observations by ~ 7 (cf Table 1). 
In this first field we detected 113 X-ray sources and in the 
second 151 (using the longest exposures in these fields). We 
detected 1 variable source in field 1 and 3 in field 2. This 



6 COMPARISON WITH SOLAR FLARES 

Flaring activity has been seen on the Sun at many wave- 
lengths and timescales. We have searched for solar flares that 
have a similar temporal profile to 'JTP 138' using the Geo- 
stationary Operational Environmental Satellites (GOES). 
The satellites are in geostationary orbit, their purpose being 
to monitor the Sun's X-ray emission continuously. 

We have identified several solar flares which have similar 
timescales to that seen in JTP 138. Figure 5 (curve a) shows 
the light curve of a 'small' solar flare which has an intensity 
several orders of magnitude lower than the largest flares in 
the current solar cycle. Like JTP 138 it shows a small in- 
crease in flux (lasting ~10 mins) before a rapid rise which 
takes around 10 mins to reach maximum flux. We have esti- 
mated its temperature using ratios of the low energy channel 
on GOES (1-8A, 1.5-12.4keV) with the high energy channel 
(0.5-4A, 3-25keV). The temperature reached a maximum of 
0.6 keV, with a luminosity of 7xl0 25 ergs s _1 . The second 
example solar flare is also shown in Figure 5 (curve b). This 
flare is brighter, and reaches a temperature of 0.8 keV, with 
a luminosity of 2xl0 26 ergs s _1 . In contrast to the small 
solar flare and JTP 138, there is no evidence for a pre-flare 
increase in flux. Both solar flares last around 1-1.5 hrs, sim- 
ilar to JTP 138. In Figure 6 we show the energy resolved 
light curve of the small solar flare. Similar to the flare seen 
in JTP 138, the solar flare reaches a maximum in hard X- 
rays earlier than in soft X-rays. 

Although the shape and duration of these solar flares 
are similar to that of JTP 138, they have lower temperatures 
and much lower luminosities. The luminosity may simply be 
a reflection of the volume of the emitting plasma. In the case 
of RS CVn systems, the flare size can extend to the size of 
the binary system (or indeed, originate in the intra-binary 
region as seen in the RS CVn system CF Tuc, Gunn et al 
1997). For reasonable system parameters and an orbital pe- 
riod of 2 days, the separation is ~ 5 x 10 11 cm. Moreover, 
the Suns rotation rate has been slowing down since it joined 
the main sequence. It therefore has lower magnetic activ- 
ity compared to stars which have recently joined the main 
sequence (as in the case of stars in NGC 2516). 

Numerous explanations for the triggering mechanism of 
solar flares have been developed over the years. The most 
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Figure 5. Two solar flares in the wavelength range (1-8A). Flare 
a) which took place on 3 Sep 1996 is a relatively small solar flare, 
which lasts for approximately 1.5 hrs. It shows prominent pre- 
flare activity. The flux has been multiplied by 10. Flare b) (12 
Feb 2001) is an average intensity solar flare. The units for flux is 
ergs/s/cm . 




2000 4000 

Time (sec) 



1.0> 
5.0> 




2000 4000 6000 8000 

Time (sec) 

Figure 6. We show the light curve of the flare a in Figure 5 
split into soft and hard X-ray bands. Like the flare seen in JTP 
138 the hard flux drops more rapidly than at softer energies and 
maximum hard flux is reached earlier than at softer energies. 



popular current theory for converting magnetic energy into 
kinetic and thermal energies is that of magnetic reconnec- 
tion (eg Priest & Forbes 1999). The site of the reconnection 
takes place high in the solar corona above magnetic loops, or 
between at least two loop systems that interact. At the re- 
connection site, the electrons are accelerated and move along 
the loops with a high velocity (v ~ 0.5c). Some of these elec- 
trons will bombard the chromosphere which is denser and 
cause this plasma to be heated. This plasma will then evap- 
orate into the loops causing the soft X-ray emission such as 
is shown in Figure 5 (cf Harra 2000 for more details). Al- 
though the details are still uncertain, observations support 
this view. 

Feldman, Laming & Doschek (1995) found a correlation 
between solar and stellar flare temperatures and emission 
measures. Shibata & Yokoyama (1999) proposed a model 
based on a magnetic reconnection model with heat conduc- 
tion and chromospheric evaporation to account for this cor- 
relation. Based on our spectral fits to JTP 138, we find an 
emission measure of 7xl0 53 cm -3 for JTP 138. For a tem- 



perature of 6x 10 7 K for the hot component, Figure 1 of Shi- 
bata & Yokoyama (1999) suggests a magnetic field strength 
of ~100 G for this source. Measuring stellar magnetic field 
strengths is difficult. However, observations of various RS 
CVn systems show evidence for magnetic field strengths of 
the order ~100-1000 G (eg Donati 1999), indicating that a 
field strength of 100 G for JTP 138 is reasonable. 



7 CONCLUSIONS 

We have performed a X-ray time variability survey of the 
open cluster NGC 2516. We have found 4 systems which 
have shown significant variability. Three of these systems 
show light curves which are characteristic of stellar flaring 
activity. Based on the flare duration and their spectral prop- 
erties we suggest that they are RS CVn systems. The fourth 
system has been identified with a AO V type star and pro- 
pose that it may have an unseen dwarf M star which shows 
flaring activity. An optical spectroscopic study of these ob- 
jects is required. Finally, this study shows that open clusters 
are excellent objects to search for stellar flaring activity and 
may be relatively common. 
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